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Abstract 

Calculations of transition energies between low-lying states of mercury atom 
are performed in the frame of combined CI/MBPT2 method. Results of 
all-electron relativistic calculations (using the Dirac-Coulomb Hamiltonian) 
are compared with experimental data and results of other four-component 
calculations. The results of the RECP calculations are compared with the 
corresponding all-electron results in order to estimate accuracy of different 
RECPs. Contributions from correlations in different shells to the calculated 
excitation energies as well as effects of basis set truncation at different orbital 
angular momenta, nuclear models, errors in gaussian approximation of the 
GRECP components are reported. Analysis of the obtained results shows that 
at least 34 external electrons of mercury atom should be correlated and the 
one-electron basis set should contain up to h angular momentum functions in 
order to reach a reliable agreement with experimental data within 200 cm -1 . 
It is concluded that correlations of the 4/ electrons can be efficiently taken 
into account for 20 electron GRECP at the generation stage. 

SHORT NAME: CI/MBPT2 calculations of mercury. 

KEYWORDS FOR INDEXING: Configuration Interaction, Many Body Perturba- 
tion Theory, correlation structure (electronic structure), excitation energies (transition ener- 
gies), Relativistic Effective Core Potential (pseudopotential), four-component calculations, 
mercury, heavy atoms. 
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I. INTRODUCTION 



During the last few years, a considerable number of papers devoted to electronic struc- 
ture calculations of heavy atoms have been published (e.g., see |]^|) that is not only due 
to experimental requirements but because they are good test systems to check or estimate 
accuracy of different approximations before using them in more expensive molecular calcu- 
lations. 

In paper H, Generalized Relativistic Effective Core Potentials (GRECPs) were tabulated 
for atoms Hg through Rn and were tested in numerical Hartree-Fock (HF) calculations by 
comparison with the all-electron Dirac-Fock (DF) and other RECP ones. However, the 
question 

of quality of the GRECPs for description of correlation effects was not clarified yet. 
Whereas paper || contains an answer on this question from theoretical point of view, the 
present paper is devoted to our correlation structure calculations for mercury atom. For this 
atom, the RECPs were generated by other groups QJ^] where the same number of electrons 
(20) was explicitly included into calculations as in the case of the GRECP [|J. Therefore, 
one 

should expect about the same computational expenses in calculations with 

all these RECPs and the comparison of their accuracy appears to be of practical interest. 

II. METHODS AND BASIS SETS 

The GRECP method was described in details in papers |§,E|||. The main distinguishing 
features of this method are presence of non-local terms with the projectors on the outercore 
(OC) pseudospinors together with the standard semi-local ones in the effective potential 
operator and generation of the effective potential components for smoothed pseudospinors 
which may have nodes f9|. 

Theory of the CI/MBPT2 method is presented in papers Q. In this method, the corre- 
lations in the valence (V) region of an atom which are the most important ones are treated 
by the Configuration Interaction (CI) method (which is able to provide excellent results for 
the small number of correlated electrons) whereas the relatively small contributions to the 
considered low excitation energies from the large number of the core- valence correlations are 
taken into account with the help of less expensive second order Many Body Perturbation 
Theory (MBPT2). 

The program (in the jj-coupling scheme) for all-electron relativistic four-component 
CI/MBPT2 calculations (using the Dirac-Coulomb Hamiltonian) was modified to make pos- 
sible two-component RECP calculations (with the nonrelativistic kinetic energy operator 
and relativistic j-dependent potentials). This program allows one to use two different basis 
sets of numerical spinors at the CI and MBPT2 calculation stages for description of the 
correlations taking into account a space separation of the core and valence regions. Basis 
functions for the present calculations were obtained from numerical SCF (DF or HF) calcu- 
lations of the corresponding spinors for positive ion states. For example, the [4,4,3,2,1] basis 
set for the 2 electron valence CI (2e-CI) was derived from the calculations for the follow- 
ing nonrelativistically averaged configurations of Hg: 5d w 6s 2 , [ScP^s 1 , [bcf^s 1 , [ScP^s 1 , 
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[bd^ds^p 1 , [5d 9 }7p\ [5d 8 }8p\ [5d 7 }9p\ [5d 9 }6d\ [5d 6 }7d\ [5d 3 }8d l , [5d 8 ]5/\ [5d 7 }6f l , 
[5d ]5g l , where the shells in the square brackets were frozen in the calculations and the ls- 
5p shells are dropped out for brevity. For convenience of comparison, a [k a , k p , k d , kf, k g , k h ] 
basis set for the case of N correlated electrons will refer here to a basis set consisting of 
n 8 sx/2, ...(n s + k s - l)si/ 2 ; n p p 1/2 , . . . [n p + k p - l)pi/ 2 ; n p p 3 / 2 , . . . (n p + k p - l)p 3 /2! n d d 3/2 , 
...(n d + k d - 1)4/25 n d d 5 / 2 , ...(n d + k d - 1)4/2! "//s/2, •••(«/ + kf - l)h/2] riff 7/2, ■ ■ ■ 
(n f + k f - l)f 7/2 ; 5g 7/2 , . . . (4 + k g )g 7/2 ; 5g 9/2 , . . . (4 + k g )g 9/2 ; Qh 9/2 , . . . (5 + k h )h 9/2 ; 6h n /2, 
. . . (5 + kh)h\i/2 spinors where n s = 5 if N > 20, n p = 5 if N > 18, n d = 5 if iV > 12, rif = 4 
if iV > 34, otherwise n s = 6, n p = 6, n d — 6, n/ = 5. Additional information on the used 
basis sets can be found on |http:/ /www.qchem.pnpi.spb.ru| . The above listed basis set has 
provided the lowest total energy in our CI calculations as compared with other tested 

similar-sized basis sets. More important advantage of this basis set is a possibility of 
one-to-one comparison between results of the all-electron 

and RECP calculations not only for large (close to full) basis set sizes 
but for the small ones as well. In fact, the calculated RECP errors depend only slightly 
on the size of such basis sets. 



III. RESULTS AND DISCUSSION 

Series of the correlation structure calculations were performed for the number of cor- 
related electrons varied from 2 to 34 and one-electron basis sets truncated at the orbital 
quantum numbers from 2 to 5. The all-electron relativistic calculations were implemented 
for two nuclear models: a point nucleus and an uniformly charged ball with 1.334 ■ 10~ 4 a.u. 
radius. The GRECP calculations were made with both the numerical GRECP components 
and their gaussian expansions from [|J. Moreover, 34 electron GRECP variant (where the 
outercore 4f,5s,5p,5d and valence 6s, 6p shells of mercury atom are explicitly treated in 
calculations) was generated and tested in these calculations. 

The results of the 2e-CI calculations (where all the possible excitations of 2 valence 
electrons were considered and the core lsi/ 2 -5ds/ 2 spinors were frozen from the SCF calcu- 
lations of the ground state with the 6s 2 configuration) are presented in tables | and [H] for 
[8, 8, 7, 6, 5, 4] and [4, 4, 3, 2, 1] basis sets. One can see from these tables that the use 

of the [4, 4, 3, 2, 1] basis set allows one to describe adequately the V-V 

correlations and that the RECP errors are rather stable in respect to the basis set size 
variation. Whereas the correlations only in the valence region are considered, the 20 electron 
GRECP is about 8 times more accurate than the 20 electron RECP of Ross et al. and 
about 30 times more accurate than the 20 electron energy-adjusted PseudoPotential (PP) J7J 
(first of all, due to the neglect of the difference between the outercore and valence potentials 
in the RECPs [|,|; see gj| for details). 

From tables III and IV in || one can see that the errors of the energy-adjusted PP, Ross et 
a/.'s RECP, and the GRECP for the excitations from the Qsf , 2 state to the 6s\ > 2 Qp\ , 2 ( J = 0) 



3 



and Qs\^ 2 6py 2 (J = 2) statesQare -667 and +347 cm" 1 , +182 and +224 cm -1 , -9 and -9 cm" 1 , 
correspondingly. They are in agreement with the errors from tables | and || here with an 
exception of the GRECP errors (because the main contribution to the GRECP errors is due 
to errors of reproducting the two-electron integrals rather than drawbacks of the effective 
potential operator). 

In tables |H|-[V|, results of both the all-electron relativistic and 20 electron GRECP calcu- 
lations are presented for different numbers of correlated electrons^ using the equivalent basis 
sets. One can see that the 2e-CI (table |H|) gives only rough description of the correlation 
structure of mercury atom, the errors in energies of excitations from the ground state are 
3000-6000 cm -1 in respect to the experimental data or the most elaborated calculations in 
table [V|. Such a large 

deviation is mainly due to neglect of the OC-V correlations involving the 5d shell. Their 
consideration (in the 12e-CI/MBPT2, table |lHj ) reduces these errors to a level of about 
1000 cm" 1 . The contributions from the correlations with the 5p and 5s shells (the 18e- 
CI/MBPT2 and 20e-CI/MBPT2, tables |TTT] and |VD to the excitation energies are 400- 
600 cm -1 and about -300 cm" 1 , correspondingly. However, these contributions are mainly 
compensated and the total contribution from the correlations with both these shells is only 
100-300 cm" 1 . The correlations with the 4/ shell (the 34e-CI/MBPT2, table |VJ) give an 
essential contribution that is 600-700 cm" 1 . Our final all-electron results for transitions 
between the first four states are within 200 cm -1 with the experimental data |T(| and the 
Relativistic Coupled Cluster (RCC) calculation results of Eliav et al. . A relatively large 
deviation for the last state is rather due to the MBPT2 approximation than due to basis set 
incompleteness or correlations in more inner shells. Our estimates show that the contribution 
from the correlations with the 4d shell to the excitation energies is of order of 100 cm -1 . 

When only the V-V correlations are considered (tables [I] and 

11), the GRECP errors in reproducting the all-electron results are within 30 cm -1 . This is 
a good confirmation to our previous estimates M that the GRECP describes the electronic 
structure in the valence region with a high accuracy. However, the consideration of the 
correlations in the outercore region (tables |T| and IV) leads to an increase of these 

errors (due to a rather large smoothing region for the nodeless outercore 5s, 5p, 5d pseu- 
dospinorsQ) up to 200 cm" 1 . These results are in a good agreement with that from table III 
in paper where a significant increase in the GRECP errors can be observed for the case of 
excitations from the 5d shell. One should expect that the contribution from the correlations 



1 Analysis for the case of the 6s\^ 2 6p\^ 2 (J = 1) and 6s\^ 2 6py 2 (J = 1) SCF states is more 
complicated because they are strongly mixed in the CI calculations. 

2 Electrons are included in the calculations in the following order: 6p, 6s; 5d; 5p; 5s; 4/; . . . 

3 We suggest that these GRECP errors can be seriously reduced due to more artificial generator 
configuration selection and smoothing procedure and with the help of corrections to the GRECP 
operator (see ||, subsection 4.5) because these errors arise from electronic structure reorganization 
in the outercore region. 
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Ill 

the outercore region to low excitation energies and GRECP errors in their reproducing 
will be reduced when passing from Hg toward Rn. 

The errors of different RECPs for the case of 20 correlated electrons are presented in 
table [TV]. The errors of the GRECP are only about 4 times less than that of Ross et al.'s 
RECP @ and about 7 

times less than that of the energy- adjusted PP J7[ for the same number of electrons 
explicitly included in calculations. This is also in agreement with the results from table III 
in paper || where the ratio between GRECP errors and errors of the other RECPs was 
reduced from one order of magnitude (when excitations only for the valence electrons were 
under consideration) to 1.5-2 times (for the case of excitations from the outercore shells 
described by nodeless pseudospinors) . In fact, additional changes in the RECP errors when 
the OC-V correlations are taken into account (tables |I| and |V|) are within 200 cm -1 for the 
used GRECP variants and within 400 cm' 1 for the other RECPs. The main reason is the 
spinor smoothing which has similar features for all these RECPs. 

The data from the MRCI calculations with the energy- adjusted PP using the CIPSO 
method for the transitions between the first four states from table 6 in paper are within 
100 cm -1 with the experimental data but the 20 electron PP does not take into account 
contributions from the correlations with the 4/ shell (which are up to 

700 cm -1 ) and the used basis set does not contain functions with h orbital momentum 
(their contributions are up to 300 cm -1 ). Therefore, these data are results from cancelation 
of a few contributions: PP errors (e.g., the 6s 1 / 2 6p{/ 2 (J = 0)-6s\/ 2 6py 2 ( J = 2) splitting is 
overestimated about 1000 cm' 1 by the energy- adjusted PP because of the features of the 
spin-orbit simulation in the LS'-based variant of the energy-adjusted scheme), a neglect of 
the correlations with the 4/ shell, a basis set incompleteness, etc. The estimation of error 



for these data within 5% appears to be correct. The data in paper [11 were obtained 
with the 12 electron RECP of Ross et al. which differs from the used here 20 electron RECP. 
One can see from tables |, U and |IV| that the errors of the 20 and 34 electron GRECPs 



are of the same order of magnitude because the main distinctions between these GRECP 
variants are 

inclusion of the 4/ electrons in calculations with the 34 electron GRECP and smoothing 
the 5/ spinors for the 20 electron GRECP. However, the 34 electron GRECP allows one to 
take into account explicitly the correlations with the 4/ shell which are important for an 
agreement with the experimental data within 200 cm' 1 . As one can see from table [V], the 
34 electron GRECP errors are within 200 cm -1 in this case. 

The basis set truncation effect at different orbital quantum numbers on energies of ex- 
citations from the ground state can be observed in tables [V] and [VI| for the 34e-CI/MBPT2 



calculations. It is clear that (spd) correlation basis set[] does not allow one to take into 
account the correlations with the 5d shell properly. This leads to the excitation energies 
(table [VTJ) rather close to that from the 2e-CI calculations. Addition of functions 



4 The [12, 12, 11, 1] basis set contains only such functions with / orbital momentum which corre- 
spond to the 4/ 5 /2 and 4/ 7 / 2 spinors. 
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with / orbital momentum is crucial for description of these correlations and leads to 
an increase in the excitation energies on 2000-4000 cm -1 . The excitation energies for the 
(spdf) basis set are within 900 cm -1 from that for the (spdfgh) basis set (tables |V| and |VT]) . 
In turn, (spdf) basis set is inadequate for description of the correlations with the 4/ shell. 
The results for this basis set are rather close to that from the 20e-CI/MBPT2 calculations 
in table |V]. Addit ion of functions with g orbital momentum is necessary for the correct 
description of these correlations and gives the contribution 300-600 cm -1 to the excitation 
energies (table |VI|) whereas addition of h functions contributes about 300 cm -1 (table |V|). 
Our test calculations showed that the contribution from functions with i orbital momentum 
is of order 50 cm -1 . 

One can see from comparison between the 6 and 7 columns in tables Q, H, and fV] that 
the errors 

due to the gaussian approximation of the GRECP components are approximately one 
order of magnitude less than the errors of the numerical GRECP. As one can see from 
tables |, |I[ |TV|, and [V], the effects of different nuclear models may be neglected for accuracy 
within 200 cm -1 . 



IV. CONCLUSIONS 

The GRECP allows one to reproduce the electronic structure in the valence and outercore 
regions essentially better than the other tested RECPs for 
the same number of explicitly treated electrons. 

At least 34 external electrons (occupying the 4/, 5s, 5p, 5d, 6s, 6p shells) of mercury atom 
should be correlated and the one-electron basis set should contain up to h angular mo- 
mentum functions in order to obtain a reliable agreement with experimental data for low 
excitation energies within 200 cm -1 whereas the errors of the gaussian approximation of the 
GRECP components and the effects of different nuclear models are negligible for this accu- 
racy. However, our test calculations show that the main contribution from the correlations 
with the 4/ shell is due to the one-electron correction from the self-energy diagrams @, 
therefore, this contribution can be taken into account for 20 electron GRECP with the help 
of the technique proposed in || (subsection 5.2). 
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TABLES 



TABLE I. All-electron and RECP correlation energies 3 , all-electron transition energies (TE) 
and absolute errors (AE) of different RECPs in their reproducting from the 2 electron valence CI 



calculations of low-lying states of Hg for the [8, 8, 7, 6, 5, 4] basis set (in cm 1 ). 

20 el. 20 el. 

Sym- Leading All-el. 34 el. 20 el. RECP energy- 

metry conf. finite point GRECP GRECP of Ross adjusted 

(Jparity) nucl. nucl. num. num. gaus. et al.° PP d 

Correlation energies 3, 

g 6s 2 1/2 -5991 -5987 -6019 -6023 -6024 -6006 -5990 

U 6s\/2 6 Pi/2 " 922 " 922 " 927 " 927 " 927 " 926 " 933 

l u 6sJ /2 6p} /2 b -6919 -6916 -6936 -6937 -6937 -6923 -6735 

2 U 6s l/2 6 ^3/2 ~ 1108 " 1108 " 1112 " im " 1116 

l u 6sJ /2 6^ /2 b 171 172 164 163 163 171 17 

TE TE AE AE AE AE AE 

g 6sj /2 

U 6s\/2 6 Pi/2 31832 31891 - 1 17 16 203 " 687 

l u §s\/2 Q Pi/2 h 33704 33764 " 3 17 16 215 " 439 

2 U 6s \/2 6 Pl/2 38073 38136 " 8 18 17 248 360 

l u §s\/2 Q P\/2 h 50108 50162 7 34 33 245 " 68 



a The correlation energies were calculated as a difference between the total energies from the above 
mentioned CI calculations and the numerical SCF calculations with the frozen core lsi/2 — 5^5/2 
spinors for the corresponding terms. 

b The 6sJ/ 2 6p]y 2 and 6s* / 2 6 P3/2 configurations are strongly mixed for these terms in the CI calcu- 
lations. 

c The RECP from Ref. §. 

d The PP from Ref. Q with the corrected V so by the factors (21 + l)/2 (M. Dolg, private commu- 
nication) . 
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TABLE II. All-electron and RECP correlation energies a , all-electron transition energies (TE) 
and absolute errors (AE) of different RECPs in their reproducting from the 2 electron valence CI 



calculations of low-lying states of Hg for the [4, 4, 3, 2, 1] basis set (in cm 1 ). 

20 el. 20 el. 

Sym- Leading All-el. 34 el. 20 el. RECP energy- 

metry conf. finite point GRECP GRECP of Ross adjusted 

(Jparity) nucl. nucl. num. num. gaus. et al. c PP d 

Correlation energies 3, 

g 6s 2 1/2 -5927 -5923 -5954 -5964 -5965 -5941 -5924 

U 6s l/2 6 Pl/2 " 909 - 910 " 914 - 915 " 915 " 913 ~ 919 

l u 6s\/2 6 Pi/2 h " 6903 - 6900 - 6920 " 6922 " 6923 " 6907 " 6718 

2 U 6sJ /2 6pi /2 -1097 -1097 -1100 -1101 -1101 -1100 -1106 

l u 6sJ /2 6^ /2 b 237 238 227 221 221 235 79 

TE TE AE AE AE AE AE 

g 6sj /2 

U 6s\/2 6 Pi/2 31780 31839 - 1 22 21 202 " 687 

l u 6sJ /2 6p} /2 b 33655 33716 -4 21 21 214 -438 

2 n 6s i/2 6 ^3/2 38020 38083 " 9 22 21 247 358 

l u 6s l/2 6 ^3/2 b 50109 50164 3 31 31 241 " 73 



a The correlation energies were calculated as a difference between the total energies from the above 
mentioned CI calculations and the numerical SCF calculations with the frozen core lsi^^c^^ 
spinors for the corresponding terms. 

b The 6sJ/ 2 6p]y 2 and ^Sy 2 6p^, 2 configurations are strongly mixed for these terms in the CI calcu- 
lations. 

c The RECP from Ref. §. 

d The PP from Ref. Q with the corrected V so by the factors (21 + l)/2 (M. Dolg, private commu- 
nication). 
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TABLE III. All-electron transition energies (TE) and absolute errors (AE) of the 20 electron 
GRECP in their reproducting from the 12 and 18 electron CI/MBPT2 calculations of low-lying 
states of Hg for the [4,4,3,2,1] CI, [11,11,11,9,8,7] and [11,12,11,9,8,7] MBPT basis sets (in 
cm -1 ). 



Sym- 




All pi 


90 el 
zu ei. 


All p1 


9n el 

zu ei. 


metry 


Leading 


point 


GRECP 


point 


GRECP 


(J parity) 


conf. 


nucl. 


gaus. 


nucl. 


gaus. 


Number of correlated electrons: 


12 


12 


18 


18 






TE 


AE 


TE 


AE 




6s i/2 














o„ 


6aJ /2 6pJ /2 


36753 


-194 


37264 


-211 


In 


6s\ /2 6p\ /2 


38526 


-178 


39049 


-192 


2 


6s\ /2 6pl /2 


43149 


-167 


43751 


-172 




6s \/2 6 Pl/2 


52919 


-30 


53325 


-24 



TABLE IV. All-electron transition energies (TE) and absolute errors (AE) of different RECPs 
in their reproducting from the 20 electron CI/MBPT2 calculations of low- lying states of Hg for 
the [4, 4, 3, 2, 1] CI and [12, 12, 11, 9, 8, 7] MBPT basis sets (in cm -1 ). 

















20 el. 


20 el. 


Sym- 


Leading 


All-el. 


34 el. 


20 el. 




RECP 


energy- 


metry 


conf. 


finite 


point 


GRECP 


GRECP 


of Ross 


adjusted 


(J parity) 




nucl. 


nucl. 


num. 


num. £ 


;aus. 


et al. a 


ppb 






TE 


TE 


AE 


AE 


AE 


AE 


AE 


% 


6*1/2 























o u 


6s\ /2 6p\ /2 


36900 


36963 


-137 


-135 


-130 


570 


-468 


1„ 


6s\ /2 6p\ /2 


38690 


38754 


-128 


-115 


-111 


581 


-233 




6s \/2 6 Pl/2 


43415 


43481 


-127 


-94 


-90 


656 


795 


lu 


6s \/2 6 Pl/2 


53004 


53062 


-22 


50 


52 


545 


303 



a The RECP from Ref. @. 

b The PP from Ref. |?J with the corrected V so by the factors (21 + l)/2 (M. Dolg, private commu- 
nication). 
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TABLE V. All-electron transition energies (TE) and absolute errors (AE) of the 34 electron 
GRECP in their reproducting from the 34 electron CI/MBPT2 calculations of low- lying states of 
Hg for the [4, 4, 3, 2, 1] CI and [12, 12, 11, 10, 8, 7] MBPT basis sets in comparison with experimental 



data and results of RCC calculations (in cm 1 ). 

RCC b 

Sym- Leading Exper. (all-el., All-el. 34 el. 

metry conf. data a finite finite point GRECP 

(Jparity) nucl.) nucl. nucl. num. 

TE TE TE TE AE 

g 6s 2 /2 

U 6s l/2 6 Pl/2 37645 37453 37569 37634 153 

l u 6s \/2 6 Pi/2 39412 39302 39361 39426 169 

2 U 6s i/2 6 ^3/2 44043 44190 44157 44224 225 

l u 6 4/2 6 P\/2 54069 55453 53553 53612 224 



a The data from Ref. fTJj]. 

b The results from Ref. || for the [27,23,21,16,10,6] basis set (in the notations of the present 
paper) and 34 correlated electrons. 



TABLE VI. All-electron transition energies (TE) from the 34 electron CI/MBPT2 calculations 
of low-lying states of Hg for different basis sets (in cm -1 ). 



Sym- 




All-el. 


All-el. 


All-el. 


metry 


Leading 


point 


point 


point 


(.Jparity) 


conf. 


nucl. 


nucl. 


nucl. 


CI basis: 




[4,4,3,2,1] 


[4,4,3,2] 


[4,4,3] 


MBPT basis: 




[12,12,11,9,8] 


[12,12,11,9] 


[12,12,11,1] 






TE 


TE 


TE 


°, 


6s 2 











0u 


Gs\ /2 6p\ /2 


37370 


37068 


32680 


lu 


6s i/2 6 Pi/2 


39151 


38798 


34530 




6s !/2 6 ^3/2 


43906 


43296 


38847 


lu 


6s l/2 6 ^3/2 


53272 


52726 


50604 
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